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Pitch  paper  pilot  is  a  computer  program  v'fich  yields  pilot 
parameters  for  a  pitch  tracking  task  and  predicts  the  pilot 
rating  of  the  aircraft  handling  qualities.  Using  Pitch  Paper 
Pilot,  optimal  SA3  gains  are  selected  for  the  fixed  form  Stability 
Augmentation  System  of  a  high  performance  aircraft  with 
structural'  bending.  This  aircraft  vss  described  in  the  Design 
Challenge  to  the  1970  Joint  Automatic  Control  Conference.  The 
final  augmented  aircraft  responses  compared  favorably  with 
desired  normal  acceleration  response  envelopes.  The  pilot 
model  in  Fitch  Paper  Pilot  is  modified  in  this  study  to  include 
pilot  lag  and  remnant  which  results  in  greater  rating  accuracy, 
although  a  few  cases  still  show  rocm  for  improvement. 
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Preface 


This  thesis  was  basically  another  in  the  series  of 
studies  that  have  been  done  on  Paper  Pilot.  Captain 
Teddy  L.  Mollis  used  the  same  procedure  in  his  thesis 
to  design  a  SAS  for  the  variable  stability  T-33.  We  used 
the  Joint  Automatic  Control  Conference  Design  Challenge 
aircraft  as  our  subject,  since  it  was  a  much  tougher 
design  problem  for  Paper  Pilot. 

Much  time  was  spent  in  the  first  half  of  this  study 
trying  u  modify  Pitch  Paper  Pilot.  Ron  Anderson,  of  the 
Air  Force  Flight  Dynamics  Laboratory,  had  discovered  that 
ratings  were  not  very  accurate  for  high  natural  frequency, 
low  damped  flight  conditions.  Following  his  suggestions, 
we  added  pilot  lag  and  remnant  to  the  pilot  model  in  Pitch 
Paper  Pilot,  in  an  attempt  to  correct  the  deficiencies. 
Although  accuracy  was  better;  it  was  still  not  good 
enough  to  eliminate  the  problem.  This  plagued  us  once 
again  in  the  Design  Challeng  1  work,  where  a  poor  design 
resulted  for  the  high  frequency,  low  damped  case.  '.Ye 
included  some  ideas  for  finally  curing  this  problem  in 
the  "Reecommendations  for  Future  Study"  section. 

We  wish  to  express  deep  appreciation  to  Major  James 
Dillow,  of  the  AFIT  Math  Department,  ror  his  constant  help 
and  advice  during  the  course  of  this  study.  'Ye  also  want 
to  thank  Mr.  Ron  Anderson  of  the  AFFDL,  for  his  valuable 
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SELECTION  OF  OPTIMAL  STABILITY  AUG:- 
FOR  A  HIGH  PERFORMANCE  AIRCRAFT 


MENTATION  SYSTEM  PARAMETERS 
USING  PITCH  PAPER  PILOT 


I.  Introduction 


Background 

The  study  of  pilot-vehicle  systens  has  been  the  subject 
of  extensive  research  in  recent  years.  Examining  aircraft 
handling  qualities  as  a  closed-loop  system,  with  the  pilot 
in  the  control  loop,  promises  to  give  more  accurate  evalua¬ 
tion  of  aircraft  flying  qualities  than  previous  evaluations 
of  the  open-loop  aircraft  alone.  Most  of  the  effort  has 
been  devoted  to  developing  an  accurate  human  pilot  model  - 
not  an  easy  task,  since  the  human  being  is  probably  the 
most  sophisticated  and  complex  control  system  in  existence. 


Pilot  opinion  rating  has  also  received  much  attention, 
notably  in  studies  by  Cooper  (3:--).  A  rating  scale  for 
pilot  opinion  of  aircraft  flying  qualities  was  developed 
from  Cooper's  work,  and  has  now  been  fairly  well  standardized 
for  flight  test  and  research.  The  Cooper  Rating  Scale, 
shown  in  Table  1  on  the  next  page,  is  a  numerical  scale 
from  1  to  10,  where  a  rating  of  1  indicates  excellent  flying 
characteristics  and  a  rating  of  10  indicates  a  totally 
unflyable  aircraft. 


Paper  Pilot.  In  1970,  Anderson,  Connors  and  Billow 
(1: --)  developed  a  promising  now  technique  of  closed- loco 
pilot-vehicle  analysis  for  pitch  tracking  tasks,  a  computer 
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program.  ca! led  "ritc.Ii  Pspi.-r  Pilot".  This  program  was  an 
exnms ion  of  prior  work  called  "Paper  Pilot"  (2:--),  a 
VTOL  hover  tracking  task  study.  Pitch  Paper  Pilot  is  a 
fully  automated  procedure  that  predicts  pilot-vehicle 
performance  and  a  pilot  opinion  for  the  pitch  tracking 
task,  based  on  the  assumption  that  a  pilot  adjusts  his 
parameters  to  minimize  his  numerical  rating  of  the  aircraft 
dynamics . 

In  1971,  Hollis  (3:--)  used  the  Paper  Pilot  program 
as  a  design  tool  to  design  an  aircraft  stability  augmenta¬ 
tion  system  (SAS)  and  optimize  its  gains  to  minimize  the 
pilot  rating  of  the  aircraft  at  given  flight  conditions. 

The  subject  aircraft  was  a  variable  stability  T-33,  an 
aircraft  not  capable  of  an  extreme  range  of  either  Mach 
number  ox*  dynamic  pressure. 

JACC  Design  Challenge.  At  the  1970  Joint  Automatic 
Control  Conference,  Rediess  and  Taylor  posed  a  flight 
control  system  design  challenge  for  a  hypothetical  air¬ 
craft  with  extreme  variations  in  flight  conditions 
The  challenge  was  to  design  a  SAS  for  this  aircraft,  whose 
flight  envelope  included  Mach  numbers  from  0.35  to  3.00 
and  dynamic  pressures  ranging  from  150  psf  to  1500  psf. 

The  aircraft  also  has  structural  flexibility  represented 
by  its  first  two  bending  modes.  This  problem  was  dealt 
with  by  Sherrard  18:--),  using  optimal  control,  V.-ard 
(9:--),  using  classical  control,  and  Paul  (5:--),  using 
self-adaptive  control  techniques. 
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Probi :2_r:  Description 

The  problem  for  this  study  was  to  design  a  stability 
augmentation  system  for  the  JACC  Design  Challenge  Aircraft. 
The  closed  loop  aircraft  response  was  compared  with 
desired  response  envelopes  as  given  in  the  Design  Challenge. 
These  envelopes  specified  limits  on  the  normal  acceleration 
and  the  derivative  of  normal  acceleration  responses.  Also, 
there  was  a  numerical  limit  on  the  rms  wind  gust  response 
to  a  specified  input. 


Approach 

The  method  used  to  design  the  SAS  was  basically  the 
Pitch  Paper  Pilot  SAS  optimization  technique  developed 
by  Hollis.  The  entire  closed  loop  system  is  shown  in 
Figure  1.  The  optimization  of  the  SAS  was  fixed-form, 
that  is,  only  the  values  of  feedback  gains  within  the 
SAS  were  optimized  with  respect  to  pilot  rating.  The 
SAS  gains  were  systematically  changed  until  a  minimum  air¬ 
craft  closed  loop  pilot  rating  was  obtained.  Those  final 
gains,  then,  were  the  optimal  SAS  gains.  This  was  done 
for  each  of  a  set  of  discrete  flight  conditions  which 
represented  the  entire  flight  envelope.  The  feedback 
vector,  x  ,  included  elevator  deflection,  S,. ,  pitch  rate,  q, 
and  normal  acceleration,  nz.  The  fixed  form  pilot  model 


Diagram  of  Closed 
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was  part  of  the  '.  itch  Paper  Pilot  routine,  which  produced 
pilot  ratings  for  the  particular  closed  loop  dynamics  in 
use.  Before  this  technique  could  he  used,  however,  certain 
^modifications  to  Paper  Pilot  were  necessary. 

As  noted  by  Anderson  (1:12),  Pitch  Paper  Pilot  did 
not  yield  accurate  pilot  ratings  for  an  aircraft  with 
high  short  period  natural  frequency  and  low  damping.  In 
this  study,  pilot  lag  and  remnant  were  added  to  the  pilot 
model  in  an  attempt  to  make  Paper  Pilot  ratinys  correlate 
with  actual  pilot  ratings  for  aircraft  with  these  short 
period  parameters.  The  level  of  remnant  was  adjusted  to 
better  match  predicted  pilot  ratings  to  actual  pilot 
ratings . 


Overview 

Since  the  first  part  of  this  investigation  concerned 
modifications  to  Pitch  Paper  Pilot  without  considering  the 
JACC  Design  Challenge,  these  modifications  were  discussed 
in  Chapter  II.  The  JACC  Design  Challenge  is  introduced 
in  Chapter  III,  and  the  digital  optimization  cf  SAS  gains 
is  summarized  in  Chapter  IV.  Analog  simulation  results 
are  the  topic  of  Chapter  V,  and  in  Chapter  VI  some  con¬ 
clusions  on  the  entire  seudy  are  drawn  and  recommendations 
offered  for  future  study. 


II . 


Prediction 


Li-  ‘_°L  "ntir.c 
I : ;  t  rod  ucti  c n 

Pilot  Rating  prediction  for  this  study  is  based  on 
a  modified  version  of  the  "Pitch  Paper  Pilot"  computer 
program,  originally  developed  by  Anderson.  Pillow  and 
Connor?  (1:--).  This  program  determines  the  pilot  model 
parameters,  closed  loop  pilot-vehicle  performance,  and 
a  Cooper  scale  pilot  opinion  rating  of  the  aircraft's 
handling  qualities  for  a  pitch  tracking  task. 

As  originally  formulated,  this  program  used  a  simple 
pilot  model,  i  ncorporatir.g  gain,  lead,  and  a  pure  time 
delay.  The  time  delay  accounted  for  both  r.euro-muscular 
lag  and  reaction  time  delay.  As  shown  in  Pigvre  1,  com¬ 
bining  this  pilot  model  with  the  SAS,  airframe,  and 
elevator  servo  models  results  in  a  complete  closed  loop 
pilot-vehicle  model,  which  is  used  to  define  a  set  of 
state  equations  used  in  the  program. 

The  procedure  for  predicting  ratings  from  this 
closed  loop  formulation  is  based  on  the  hypothesis  that 
a  human  pilot  will  adjust  his  parameters  (the  equivalent 
of  the  pilot  model  gain,  lead  and  time  delay)  to  allow 
him  to  give  the  lowest  (best)  rating  to  the  aircraft 
that  he  can  at  the  particular  flight  condition.  Thus, 
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i  —  3 

by  minimising  seme  properly  chosen  rating  expression  with 
respect  to  pilot  parameters ,  the  program  can  be  expected 
to  produce  ratings  veil  correlated  with  actual  pilot 
ratings.  In  the  original  program,  an  empirically  developed 
rating  function,  PR,  is  minimized  with  respect  to  Kp,  the 
gain  of  the  pilot  in  the  pitch  loop,  and  T^,  the  pilot 
lead  in  the  pitch  loop  (the  time  delay  is  considered  a 
constant)  .  Conjugant  gradient  and  Uewton-Raphscn  search 
procedures,  as  shown  in  Hollis  (3:17-23),  were  used  to 
find  the  minimum. 


Rating  Expression.  Anderson,  Dillow  and  Connors 
(1:4-5)  developed  a  rating  expression  for  use  with  the 
above  procedure.  This  expression  is  of  the  form 


PR  =  Wj_  +  Vi  2  +  1.0  (1) 

where 


and 


0.1  <7C _ 

0 . 9  7  4  oc  -  a  t 


®t<  0.974ae  (2) 


W 
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-2.5  Tr 
2.5  Tl 
3.25 


Tl  <  0 

OSTl-1*3  (3) 

TL>1.3 


is  a  measure  of  closed  loop  performance  where  vc"  is 
the  commanded  input  variance  and  o~  is  the  pitch  error 
variance.  W2  is  a  measure  of  pilot  workload,  and  is  a 
function  of  T^,  the  pilot  lead. 


n 


As  orginally 
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Short.com i  ru: :■  of  tho  Orl"  Inal  JPrograrn . 
formula tod,  Pitch  Paper  Pilot  gave  accurate  ratings  for 
most  flight  conditio. .-3.  uowover,  when  the  program  was 
applied  to  flight  conditions  with  high  short  period 
frequency  and  low  damping,  the  predicted  ratings  for  the 
flight  conditions  v:ere  significantly  better  (lower)  than 
actual  ratings  (1:12). 

Since  the  Paper  Pilot  concept  does  work  very  well  for 
most  pitch  cases  and  for  all  flight  conditions  when  applied 
to  other  tasks,  such  as  roll  (5:--)  and  hever  tracking 
(2:--),  some  modification  to  Pitch  Paper  Pilot  should 
correct  this  fault. 


Remnant  Studio s 

One  possible  approach  to  improving  the  rating  scheme  i 
to  include  pilot  remnant  in  the  pilot  model.  Remnant  is 
that  portion  of  a  pilot's  output  that  is  not  linearly 
correlated  with  the  systera  input  (4:--).  Pure  noise  injec¬ 
tion  by  the  pilot,  non-linear  operations,  and  non-steadv 
pilot,  behavior  are  considered  to  be  the  possible  sources 
of  remnant. 

It  was  hypothesized  that  the  addition  of  remnant  to 
tbe  pilot  model  would  tend  to  produce  higher  ratings  from 
Pitch  Paper  Pilot  for  all  flight  conditions,  taut  that 
this  effect  would  be  greater  for  conditions  with  low 
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damping  and  high  frequency.  This  would  allow  the  critical 
(high  frequency,  lew  camping)  predicted  ratings  to  be 
raised  to  more  realistic  values  without  pushing  other 
cases  too  far  above  their  original  values. 

Considering  remnant  as  pilot  induced  noise  with  a 
relatively  high  bandwidth,  it  is  easy  to  see  that  this 
pilot  noise  could  excite  an  undesired  response  in  the  air¬ 
craft.  An  aircraft  with  a  high  short  period  natural 
frequency  should  filter  less  of  this  remnant  effect  out, 
resulting  in  greater  tracking  error.  If  this  undesired 
response  were  also  lightly  damped,  it  would  be  even  more 
noticeable  to  the  pilot.  Thus,  for  the  critical  case, 
tracking  error  and  pilot  workload  would  increase  more  than 
for  other  cases.  For  a  pilot  model  with  remnant,  therefore 
predicted  ratings  would  be  expected  to  increase  by  a 
greater  amount  for  the  critical  case  than  for  other  cases. 

New  Pilot  Model.  To  test  this  hypothesis,  a  new 
pilot  model  was  used  '.n  the  computer  program  which  was  a 
higher  order  model  including  both  remnant  and  a  lag  term. 
The  new  pilot  model  is  shown  in  Figure  2.  Note  that  the 
remnant  is  modeled  as  a  zero-mean  gaussian  white  noise 
put  through  a  first  order  filter.  To  model  high  frequency 
remnant,  the  cut-off  freqency,  u>r ,  was  chosen  as  20  radians 
per  second.  Tj  is  the  lag  time  constant,  which  was  chosen 
as  a  constant  of  0.1  second:,  to  represent  the  no uro-muscul a 
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lag.  t  in  the  time  delay  term  was  reduced  to  0.32 
seconds,  representing  the  higher  order  neuro-motor  delay 
and  the  pure  time  delay.  The  program  still  optimizes 
with  respect  to  and  only. 

A  set  of  system  equations  for  this  pilot  model, 
including  the  actuator,  are  developed  in  Appendix  B.  By 
combining  these  equations  with  the  longitudinal  airframe 
and  control  system  equations  of  Appendix  A  and  reducing 
them  to  state  variable  form,  an  accurate  mathematical  node, 
of  the  closed-loop  system  was  substituted  into  the  proper 
sub-routines  of  the  original  Pitch  Paper  Pilot.  This 
modified  program  was  then  used  to  test  the  effect  of  the 
remnant  on  rating  predictions. 

Results  of  Remnant  Studies 

The  modified  Pitch  Paper  Pilot  program  was  tested  on 

The  data  frc;v.  a 


cases  for  which  real  ratings  were  known. 
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series  of  Navy  tests  with  a  variable  stability  F-8  (1:10-12) 
were  used  as  test  cases.  Pertinent  data  for  these  two 
sets  of  cases ,  one  set  for  landing  approach  and  one  for 
low  altitude  high  speed  tests.-  are  shown  in  Tables  II  and 
III  on  the  following  pages.  Also  included  in  the  Tables 
are  the  actual  average  pilot  rating.  As  can  be  seen 
from  Table  III,  the  agreement  between  actual  and  predicted 
ratings  is  particularly  peer  for  the  low  altitude  high 
speed  cases.  These  cases  are  mainly  of  the  high  frequency, 
low  damping  type.  For  the  other  cases,  as  shown  in  Table 
II,  Paper  Pilot  did  give  good  results. 

The  Effect  of  Remnant.  Using  the  new  Pitch  Paper 
Pilot,  u  series  cf  predicted  pilot  ratings  were  computed 
for  each  flight  condition  using  a  variety  of  different 
levels  of  remnant,  starting,  in  each  case,  with  a  zero 
remnant.  A  commanded  input  with  an  rms  value  of  0.5 
radians  was  used  in  all  cases.  With  this  input  and  the 
remnant  specified,  stick  input  could  be  calculated  for 
each  flight  condition  and  remnant  level. 

Remnant  only  begins  to  affect  ratings  as  it  becomes 
significant  in  comparison  to  the  stick  input  required 
for  the  particular  flight  condition.  For  this  reason, 
the  set  of  ratings  for  each  flight  condition  were  plotted 
with  respect  to  the  remnant  to  stick  ratio  (  <rr2/  af  )  . 
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Average  pilot  rating  for  the  test  aircraft 
at  trie  given  flight  condition.  Individual 
pilot's  ratings  varied  as  rush  as  thre 
ratings  for  some  flight  conditions. 

Ratings  calculated  Iron  the  original  Pitch 
Paper  Pilot  program. 
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Plots  of  remnant  to  stick  ratio  versus  predicted 
rating  are  very  siniiar  for  all  cases;  the  predicted 
rating  rises  by  a  op  vox  i  mace  ly  the  same  amount  for  a  given 
remnant  to  stick  ratio,  no  matter  what  the  flight  condition. 
Four  typical  cases  are  shown  in  Figures  3  and  4.  They 
clearly  show  this  similarity  of  curves  for  different 
flight  conditions. 

This  similarity  in  remnant  effects  makes  it  impossible 
to  pick  a  particular  remnant  to  stick  ratio  that  would  in¬ 
crease  the  ratings  for  the  critical  case  (high  frequency, 
low  damping)  without  increasing  the  ratings  in  the  other 
cases  by  about  the  same  amount.  Thus,  the  critical  cases 
could  not  be  adjusted  to  more  closely  correlate  with  real 
pilots'  ratings  without  having  a  corresponding  effect  of 
increasing  the  other  ratings  to  values  that  were  signifi¬ 
cantly  higher  than  their  actual  ratings. 

Modified  Rating  Expression 

Including  remnant  in  the  pilot  model  failed  to  provide 
the  anticipated  improvement  in  Pitch  Paper  Pilot,  so  modi¬ 
fication  of  the  rating  expression  was  tried.  Small  adjust¬ 
ments  were  made  to  the  constants  in  the  pilot  rating 
expression  in  an  attempt  to  correct  the  critical  ratings. 
This  modification  also  failed  to  solve  the  problem.  Each 
flight  condition  was  similarly  affected  by  equal  changes 
in  constants,  as  with  the  remnant  ef::ects. 
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Final  .Mcdif icnticn  to  Paper  Pilot 

Obviously,  the  changes  made  in  the  original  program 
were  not  very  successful  in  eliminating  Paper  Pilot's  main 
fault,  inaccurate  ratings  for  high  frequency,  low  damping 
conditions.  However,  the  modified  pilot  model,  with 
remnant  and  lag,  w~s  incorporated  in  the  closed-loop  model 
because  this  model  is  a  closer  approximation  to  the  human 
pilot  and  the  use  of  a  non-zero  remnant  did  slightly 
improve  the  predicted  rating  for  the  critical  flight 
conditions . 

A  remnant  to  stick  ratio  of  0.6  was  selected  for  the 
revised  mode].  This  ratio  is  in  a  region  that  does  not 
cause  extremely  large  changes  in  rating  for  small  changes 
in  the  remnant  to  stick  ratio.  The  scatter  diagram  of 
actual  versus  predicted  pilot  rating  for  this  value  of 
remnant  to  stick  ratio,  as  shown  in  Figure  5,  does  show 
a  fair  cprrelation  of  actual  to  predicted  pilot  rating. 

The  0.6  remnant  to  stick  ratio  did  give  the  best  scatter 
diagram  of  the  ratios  tested. 

Calculation  of  Remnant.  The  remnant  required  to 
generate  the  proper  remnant  to  stick  ratio  is  calculated 
by  a  simple  procedure  incorporated  in  the  computer  program. 
Mean-square  remnant  value  is  determined  by  manipulating  the 
following  equation 

*/=  Kl"i2  4  0,2 

-1  p 


(4) 
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Predicted  Pilot  Rating 

Figure  5.  Scatter  Diagram  for  oj?A2  of  0,0, 

where  «r.2  is  the  mean-square  input.  First,  k-^  is  found 
by  setting  the  remnant  at  zero  ( <rr  =  0  )  so  that 

*i  - 

i 

then  k2  is  found  by  setting  the  remnant  at  unity  and  the 
input  to  zero  (  «r  =  0)  so  that 


ko  = 


Now,  for. a  specified  remnant  to  stick  ratio,  k-j , 


«7t2(l/k3)=  k^aj2  +  k2ar 


r  "  0/k3)-k2 

This  routine  was  added  to  Pitch  Paper  Pilot  and  the 
modification  was  used  in  the  out imizat ion  routine. 
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Though  the  modifications  were  not  totally  successful, 
the  resultant  Pitch  Paper  Pilot,  including  remnant,  was  an 
^improvement,  and -assumed  to  be  a  good  tool  for  SAS  design. 
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in.  joi nt  Autoir.av.ic  Control  Ccnfcrc^o  Dcp igr._Ch^I_lonn£. 

The  JACC  Design  Challenge  (7:--)  specified  aircraft 
dynamics,  servoactuator  dynamics,  onboard  measurement 
systems,  structural  response  dynamics,  a  mission  profile, 
a  gust  environment,  and  desired  closed  loop  responses  for 
the  final  augmented  aircraft  design.  The  aircraft  was  to 
be  piloted  by  a  human  pilot  through  the  mission  profile. 
Sixteen  discrete  points  in  the  mission  profile  were  selected 
as  representative  of  the  entire  flight  envelope.  These 
flight  conditions  are  listed  in  Tabic  IV.  The  problem  was 
to  design  a  SAS  for  this  aircraft  which  brought  the  air¬ 
craft  response  in  all  flight  conditions  to  within  the 
specified  envelopes. 

Aircraft  Dynamics 


Only  the  longitudinal  mode  of  the  aircraft  was  consi-  ’ 
dered.  The  phugoid  mode  was  assumed  to  be  negligible.  The 
longitudinal  short  period  rigid  body  dynamics  of  the  basic 
aircraft  are  given  by  the  following  set  of  linearized 
differential  equations 
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trajectory  at  fifteen  second  intervals,  with 
Flight  Condition  1  at  zero  tirrn. 
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where  a  ii  angle-of-attack ,  q  is  pitch  rate,  $ ‘  is 
elevator  deflection,  w0  is  a  wind  gust#  U0  is  nominal 
airspeed,  a  is  the  acceleration  of  gravity  and  the  coef¬ 
ficients  of  th-.-sc  variab’-'  a  are  the  aircraft  stability 
derivatives.  These  equations  are  transformed  into  a 
state  variable  form  used  in  this  study  in  Appendix  A. 

in  the  Design  Challenge,  the  stability  derivatives 
were  expressed  statistically  as  random  variables  T  Lh  a 
normal  distribution.  However,  to  narrow  the  design 
problem,  only  the  mean  values  of  the  stability  derivatives 
were  used  in  this  study.  They  are  listed  in  Table  V. 

The  original  statistical  description  of  these  coefficients 
can  be  found  in  the  JACC  Design  Challenge  (7:13). 


Measurement  System  and  Structural  Rest-: 


Dynamics 


The  effect  of  the  structural  response  on  the  feedback 
system  is  a  function  of  the  sensor  locations.  Only  the 
first  two  bending  modes  were  considered,  modeled  by  the 
differential  equations 


+  2$i“i  S  +  =  *st 

«2  +  25 


32  2  2 


+  = 


(12) 
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where  *  =  ~40oZ^  and  c, 
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Challenge  presented  the  bending  mode  parameters  in  a 
statistical  sense,  only  the  mean  values  were  used  in 
this  study.  These  are 

=  30  rad/sec  5=  0.01  (14) 

u.=  50  rad/sec  $  =  0.01 
2 

The  measurement  sensors  vere  assumed  to  be  ideal, 
however,  the  quantities  measured  included  the  bending 
mode  response  and  inc  ependent  noise  in  addition  to  the 
rigid  body  response.  The  following  measurements  were 
used  for  the  selected  SAS  configuration 


C  i  *  -  V, 


A  eJ^2  +  qn 


nz.  =  nz  +  |  (Lzq  +  0n1  +  0n2^2  ]  +  nzn  (16) 

These  are  the  indicated  pitch  rate  and  indicated  normal 

acceleration  equations.  A-  and  Ag  are  the  respective 

bending  mode  slopes  at  the  sensor  locations  relative  to 

the  generalized  displacements,!^  and  u2,  while  0n^  and 

0n^  are  the  relative  displacements  at  the  accelerometers 

due  to  the  bending  modes.  Lz  is  the  distance  from  the 

aircraft  center  of  gravity  to  the  accelerometers.  The 

terms  n_  and  q  are  zero  mean  craussian  white  noise. 

‘n  n 

The  values  of  the  parameters  used  are 


=  -0.02a  rad/ft. 

Af  =  0.04  0  rad/ft. 
(J2 
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10  ft. 

-0.15 
0.45 

(O.OOOS)^  (rad/sec)^ 

(0.01)2  (g)2 

Se  rvoactuator 

The  control  servoactuator  is  described  by  a  first 
order  lag  with  an  output  rate  limit.  The  equation  for 
the  servoactuator  is 

Se  =  -20  St  +  20  S'  (17) 

with  a  rate  limitation  of 

I  5e  I  -  0,5  rad/sec  (18) 

where  St  is  a  commanded  input  to  the  servoactuator. 

Design  Specifications 

The  design  specifications  from  the  Design  Challenge 
are  the  response  envelopes  shown  in  Figure  6.  These  are 
the  envelopes  superimposed  on  the  plots  of  the  aircraft 
responses  shown  in  Appendix  C.  The  response  envelopes  for 
normal  acceleration  ana  its  derivative  are  normalized 
in  amplitude  by  the  steady  state  normal  acceleration 


L  = 
z 

= 


E(q„|2  = 

-<"z  )2  - 
n 
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r-- 


and  in  tine  by  a  factor  t  which  is  a  discrete  function 
of  dynamic  pressure 


t  =  3.C  sec  150  S  q  <  500  psf 

r  ---•  1.0  see  500  5  q  5  1500  psf 

In  the  absence  of  a  pilot  input,  the  aircraft  respon 
to  a  specified  vertical  wind  gust  must  be  less  than  0.1  c 
rras.  The  wind  gust,  w*g ,  is  assumed  to  have  a  zero-mean 
gaussian  distribution,  and  a  power  spectral  density  given 
by 


G  ( to )  =  - 

W  71 


(19) 


where  a  =  5.0  fps  and  u>c  =  U  /1000  rad/sec. 


Finally,  the  Design  Challenge  specified  restrictions 
on  stick  force,  but  this  requirement  was  not  considered 
in  this  study.  It  vas  assumed  that  stick  force  could  be 
artificially  induced,  and  would  not  limit  the  SAS  design 
involved  in  this  study. 
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IV.  SAT'  Parameter  Optimisation 

Designing  a  stability  augmentation  system  for  the 
JACC  Design  Challenge  .aircraft  was  the  main  project  for 
this  study.  This  was  done  by  changing  the  feedback  gains 
for  a  specified  fixed  form  SAS  until  the  closed  loop 
system,  including  the  pilot,  actuator  and  SAS,  as  well 
as  the  aircraft  dynamics,  had  the  lowest  pilot  rating. 

This  optimum  pilot  rating  was  determined  by  Pitch  Paper 
Pilot. 

The  design  was  restricted  to  the  longitudinal  axis, 
for  a  pitch  tracking  task.  A  detailed  diagram  of  the  closed 
loop  system  model  is  shown  on  the  following  page. 

Figure  7.  Note  that  the  stability  augmentation  system 
used  the  three  parameters;  elevator  deflection, 
pitch  rate,  q,  and  normal  acceleration,  nz.  As  discussed 
in  Chapter  III,  measurements  of  nz  and  q  are  subject  to 
disturbance  by  the  aircraft  structural  bending.  For  this 
reason,  these  measurements  are  filtered  in  the  feedback 
loop,  to  remove  some  of  this  disturbance. 

For  this  aircraft,  the  break  frequency  of  the  filter 
was  chosen  as  8.  This  was  the  value  arrived  at  by  Sherrard, 
and  represented  a  good  choice  of  maximum  attenuation  of 
body  bending  effects  with  minimum  attenuation  of  the  rigid 
body  dynamics  (8:28). 


Detailed  Block  Diagram  of  Closed  Loop  System 


rr 


The  pilot  model  for  the  closed  loop  system  is  the 
new  pilot  model  developed  for  Pittn  apor  Pilot  in  Chapter 
IT.  The  servo  actuator  dynamics  are  represented  by  a 
simple  first  order  lag. 

Brief  Program.  Description 

Pitch  Paper  Pilot  uses  the  closed  loop  system  model 
to  determine  the  aircraft's  predicted  pilot  rating.  A 
pattern  search  technique,  developed  by  Zuckerman  (10: — ), 
uses  Pitch  Paper  Pilot  to  find  the  SAS  gains  that  minimize 
predicted  pilot  rating  within  the  actuator  rate  constraints 
given  in  Chapter  III. 

The  procedure  is  as  follows:  SAS  gains  are  changed 
systematically  and  Pitch  Paper  Pilot  computes  a  predicted 
pilot  rating  for  each  set  of  gains.  The  program  proceeds 
until  a  set  of  gains  is  found  that  give  a  minimum  predicted 
rating,  without  allowing  the  quantity,  ”  Sc  ,  to  exceed  0.5 
radians  per  second,  is  a  three  sigma  measure  of 

elevator  rate,  so  there  is  only  a  very  low  probability 
that  the  rate  constraints  would  be  violated  if  is 
kept  within  limits. 

This  SAS  optimization  technique  is  basically  the 
same  as  that  used  by  ilollis  in  his  study  of  the  variable 
stability  T-33  (3:25-31). 
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New  A  **?•-. r lx.  Pitch  Paper  Pilot  requires  the  closed 
lc op  syster.  tz  be  modeled  in  a  state  vector  equation  form. 
This  is  written  as  •  ~  ...  .  .;.T.  ......  ...... 

z  =  Az  +  v  (20) 

where  2  is  the  state  vector,  v  is  a  noise  vector  and  A 
is  a  matrix  cf  coefficients  called  the  A  matrix. 

In  examining  the  Design  Challenge  aircraft,  two 
closed  loop  models  were  used,  one  of  ninth  order  without 
the  two  second  order  structural  bending  modes  and  one  of 
thirteenth  order  including  the  bending  modes.  Thw  state 
equations  for  these  two  models  are  developed  in  detail  in 
Appendices  A  and  B. 


The  fortr.  of  the  state  vector,  z,  for  ninth  or  thirteenth 
order  is  as  shown  below. 


z 


9 

q 

P 

w 

*1 

*2 

ec 

_ 


or 


(21) 
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The  two  vectors  are  similar,  except  that  the  thirteenth 

order  includes  four  extra  states  to  represent  the  bending 

modes.  9  is  the'  pitch.  ::  is  pitch  rate.  p  is  a 

feedback  state  ccr:iiniT,g  vertical  acceleration  and  pitch 

rate  feedback  (  fi  =  K  qf  +  K  n  ) .  w  is  a  state 

q  nz  zf 

formed  from  angle  of  attack  and  true  airspeed  (w  =  U0o) . 
Sc  is  the  actuator  (elevator)  deflection.  y,  and  y ^  are 
states  formed  from  the  lag-lead  and  time  delay  portions 
of  the  pilot  model,  as  detailed  in  Appendix  B.  ©c  is 
the  pitch  command.  <5r  is  the  remnant.  Uj  and 
represent  the  longitudinal  bending  terms  for  the  JACC 
aircraft  and  r  and  s  are  their  respective  derivatives. 


The  A  matrices  for  these  models,  which  are  necessary 
for  Pitch  Paper  Pi-.-'t-  are  shown  on  the  following  pages. 
Figure  8  shows  the  ninth  order  A  matrix  and  Figure  9 
shows  the  thirteenth  order  A  matrix. 


0£  timization  Results 

The  SAS  gains  were  first  optimized  for  the  ninth 
order  model,  with  no  structural  bending  considered  for 
fifteen  flight  conditions.  Table  VI  shows  the  results. 

The  table  shows  the  minimum  predicted  rating,  the  aug¬ 
mentation  system  gains  that  produced  it,  and  the  calculated 
three-sigma  value  of  actuator  rate  required.  Note  that 


Figure  8.  “a"  Matrix  without  Bending  Modes 
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A  31."  V 


Fli-ht 

Condition 


-0.24 

-0.17 

-0.09 

-0.04 

-0.04 

-0.04 

-0.07 

-o.on 

-0.06 

-0.09 

-0.11 

-0.15 

-0.17 

-0.20 

-0.04 


»  *  <  '***  j  \  *23  *  m  r 

,  i  ;  .  \  O  |  .  .A  ■»*.<.»  4  .  f  V  -  . 

irl  7IK7H  ORDAR  MCOi 


Kr. 

(r  ar./a) 


-0.0014 

-0.0002 

-0.0001 

0.0001 

0.0000 

0.0001 

0.0002 

0.0002 

0 

0.0000 
0.0004 
0.0006 
0.0007 
0.0004 
C  .001  7 


-0.63 

-0.20 

0.09 

-0.3B 

-0.48 

-0.34 

-C.54 

-0.67 

-0.49 

-0.59 

-0.74 

-C.77 

-0.36 

-0.63 

-0.17 


Actuator 

Rato 

(rad/soc ) 

Ratir^ 

0.493 

2.900 

0.493 

2.597 

0.499 

2.510 

0.499 

2.447 

0.499 

2.45? 

0.499 

2.325 

0.495 

2.184 

0.499 

2.224 

0.497 

3.263 

0.487 

2.916 

0.498 

2.415 

0.499 

2,620 

0.499 

3.972 

0.499 

2.808 

0.499 

3.474 
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FOUND  IN  THE  NINTH  ORDER  OrTI.VIZATICN. 


Flight 

Ccn<3  i  t  ion 

Actuator 

Rate 

(rad /see) 

F  r  e  "  i  c  -  -o 
Rati ag 

1 

.502 

2.900 

2 

.501 

2.590 

3 

.500 

2.51 1 

4 

.500 

2.448 

5 

.501 

2.450 

6 

.500 

2.328 

7 

.499 

2.184 

8 

.502 

2.222 

9 

.501 

3.262 

10 

.489 

2.925 

11 

.500 

2.412 

12 

.503 

2.61  5 

13 

.501 

3.975 

15 

.502 

2.803 

16 

.502 

3.475 

t i  i,  *1  , ,  .  -»  'i .  %  (  :  1  i .  •  «  -  r  »<"  ,'i  I,  it 

*• 

r  *  ,-m  i  ♦  *  i ,  i  *.  1*  ■  4* <  •  h  i  f*  li  i  *  #  *  r  «11  flight 

ev*  iit  ic.  »  . 

West,  the  thittt-enih  o»<>r  systeii  was  run  without 
f  A  3  opT'’l'tMti  *i  ,  ini^T  th*  cpt l~  » 1  feedbat’*  gain*  from 
t 1  •  Mn*h  tfi^f  i  *>■ »  •  1  ♦  *  .  A*  cun  i  «  from  the  results, 

taiul***)  in  *  hi*  '  II,  both  ;t*li  ♦  I  tat.;*  ar.2  set :.‘:r 
rates  wi*  clos.iy  (•)«.{  »[»!  le  f'  r  ooth  r  o  i*  1  s  ,  thou  gh 
actuator  rat*  limitation*  were  exceeded  for  iom  cases. 
Apparently,  the  bervlin;  modes  have  little  effect  on  th* 
design  in  these  coses, 

>a  *  further  tost,  selected  flight  c.n-Htions  with 
varied  Mach  n>imber  end  dynamic  pressvres,  were  fully 
Optimized  with  the  thirteenth  order  rod*  l  to  see  how 
th*  optimal  SAS  gains  were  changed.  Table  VIII  shows 
th*  results  "of  this  test.  Again,  only  negligible  variation 
from  th*  ninth  order  results  was  observed. 


Table  VIII.  Optimum  Fsra-»teri  for  Selected 
flight  Conditions  with  the  Thirteenth  Order  Itode  1 . 


riight 

Condition 


*n 

(ado) 


(rad?}) 


-0.24  I  -0.0014 
0.5301 


0 . 0  f,  o  'j 


Act  ua  tor 
Bate 

(rad/sec) 

Rating 

0.409 

2.901 
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For  the  J.'CC 'Design  Challenge,  addition  of  the 
bending  nodes  made  little  difference.  Ninth  order  results 
'could  be  used  for  aL  huditiors .  This  has  an  extra 
advantage,  since  much  loss  conputer  time  was  required  for 
ninth  order  optimization. 


0 


C 


Gain  Scheduling 

In  his  study  of  a  variable  stability  T~33,  Dollis 
found  a  simple  relationship  between  l/c^  and  the  SAS 
feedback  gains  (3:33-38).  This  allowed  the  gains  to  be 
scheduled  for  the  flight  conditions,  in  relation  to 
dynamic  pressure. 

The  JACC  Design  Challenge  optimization  results 
showed  a  similar  correlation  between  l/o_  and  pitch  rate 
feedback  gain.  This  is  shown  in  Figure  10.  Note  that  the 
gains  are  all  grouped  along  a  curve,  with  the  exception  of 
Flight  Condition  16.  No  reason  was  found  for  this  single 
variation . 

No  simple  dependency  could  be  found  for  elevator 
deflection  feedback  gain,  .  Figure  11  shows  a  graph 
of  K^.  versus  1/^. 

Normal  acceleration  feedback  gain  was  so  close  to 
zero  in  all  cases  that  it  would  probably  be  considered 
zero  for  scheduling  purposes. 


••  n 


-  --  •  ‘Addition  of  Pre-Filter 

Analog  simulation  of  the  Design  Challenge  aircraft 
and  S'.S  ,  vs  presented  in  Chapter  7,  vovealed  that  a  pre¬ 
filter  on  the  pilot's  input  to  the  aircraft  was  needed  to 
meet  the  design  criteria.  Without  the  pre-filter,  an 
undesirable  negative  spike  was  present  in  the  initial 
normal  acceleration  and  acceleration  rate  which  exceeded 
the  envelopes  defined  in  Chapter  III,  The  pre- filter 
significantly  improved  this  response  for  most  flight 
conditions. 

Adding  the  pre-filter  produced  a  new  system  model 
of  fourteenth  order.  Since  the  pre-filter  was  not  in 
the  original  system  model,  it  was  reasonable  to  expect 
that  the  new  model  would  have  seme  effect  on  the  predicted 
ratings  from  Paper  Pilot. 

The  pre-filter  is  on  the  ,fp  signal  from  the  pilot 
model  to  the  actuator  (see  Figure  7).  Thus,  the  model 
is  modified  as  shown  in  Figure  12. 

The  break  frequency,  u>p  ,  used  here  was  5  rad/sec. 

The  new  state  equations  for  this  fourteenth  order  mode 
are  fully  developed  in  Appendix  D. 

Table  IX  shows  the  results  from  Paper  Pilot  studies 
with  the-  new  state  equations.  The  feedback  gains  from  the 


ninth  order  SAS  optimization  were  used  to  find  new  predicted 
pilot  ratings  and  actuator  rates.  As  shown  in  Table  IX, 
the  addition  of  the  pre-filt.er  produced  a  substantial 
change  in  botu  predicted  rating  and  actuator  rate.  Note 
that  all  ratines  were  lowered  by  about  the  same  amount. 

With  both  the  rating  and  actuator  rate  lowered,  it 
was  expected  that  a  SAS  gain  optimization  for  any  flight 
condition  would  further  lower  predicted  pilot  ratings  while 
raising  the  actuator  rate  close  to  the  limit.  Flight 
condition  1  was  checked  to  test  this  hypothesis.  As  shown 
in  Table  X  below,  the  expected  result  occurred. 


Although  sore  studies  were  done  with  tho  system 
including  the  pro- filter,  a  great  deal  more  work  could 
-be  dona  in  this  area. 

Conclusions 

The  SAS  optimization,  procedure,  including  Pitch 
Paper  Pilot,  worked  well  for  this  system.  In  general,  tho 
procedure  works  well  for  systems  of  a  reasonable  order  and 
with  few  feedback  states.  Increasing  the  order  of  the 
system  or  the  number  of  feedback  parameters  requires 
causes  a  large  increase  in  computer  time.  For  instance, 
increasing  the  order  from  ninth  to  thirteenth  order 
approximately  doubled  the  required  computer  time.  The 
Air  Force  Flight  Dynamics  Laboratory  has  improved  some 
of  the  search  routines  within  Paper  Pilot  to  reduce  time 
requirements,  which  should  help  to  alleviate  this  problem. 
As  stated  earlier,  the  bending  modes  did  not  significantly 
influence  the  SAS  design  in  this  case. 
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V.  Simulation 


Procedure  -  Commanded  Se  Response 


The  Design  Challenge  aircraft  and  the  SAS  were 
simulated  on  an  EAI  Pace  analog  computer.  For  comparison 
purposes,  and  to  check  the  computer  patching,  the  system 
was  first  run  open  loop,  that  is,  without  the  stability 
augmentation  system.  Then,  the  feedback  loops  were  closed. 
The  input  for  both  cases  was  a  commanded  step  elevator 
deflection.  The  plots  are  presented  in  Appendix  C. 


Throughout  the  simulation,  only  the  mean  values  of 
the  stability  derivatives  and  flexure  mode  parameters 
were  used.  The  additive  white  noise  terms  in  the  equations 
of  the  states  perturbed  by  the  flexure  were  dropped.  It 
was  assumed  that  the  analog  computer  generated  enough  inner 
noise  to  compensate  for  the  elimination  of  these  terms. 


Flight  Conditions  1,  3,  5,  9  and  13  are  presented  here 
as  a  cross-section  of  the  possible  combinations  of  Mach 
number  and  dynamic  pressure.  Table  XI  gives  a  comparison 
of  these  values  for  the  five  flight  conditions. 


JACC  Design  Challenge  Response Envelopes 


Before  examining  the  closed  loop  aircraft  responses, 
several  points  should  be  made  about  the  JACC  Design  Challenge 
normal  acceleration  response  envelopes.  Their  shapes  appear 
to  be  contrived  limits  on  more  conventional  parameters  sucn 
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Table  XI.  Comparison  of  Flight  Conditions 


Flight 

Condi  tier. 

Dynamic 

Pressure 

Mach 

Number 

1 

low 

low 

3 

medium 

medium 

5 

high 

medium 

9 

high 

high 

13 

low 

high 

as  rise  time,  settling  time,  peak  overshoot  and  under¬ 
shoot.  It  is  difficult  to  say  whether  or  not  these 
envelopes  actually  represent  good  handling  qualities  in 
the  opinion  of  a  pilot.  Consequently,  the  Design  Challenge 
envelopes  may  not  be  good  criteria  for  evaluating  the 
response  of  an  aircraft  with  a  SAS  design  based  on  pilot 
rating . 

In  addition,  the  normalizing  factor,  t  ,  was  a 
discrete  function  of  dynamic  pressure.  Depending  on  how 
it  was  arrived  at,  this  non-linearity  could  over¬ 
restrict  responses  for  eitner  the  extreme  values  of  g_ 
or  for  the  middle  values.  In  any  case,  r  as  a  continuous 
function  of  2.  would  have  perhaps  been  a  better  normalizing 
constant . 

Regardless,  the  Design  Challenge  response  envelopes 
were  used  as  the  criteria  to  judge  the  normal  accelera- 
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tion  response  of  the  augmented  aircraft.  They  did  pro¬ 
vide  a  quick  "eyeball"  assessment  of  the  final  closed 
lo op  response. 


Normal  Acceleration  Response 

All  five  flight  conditions  had  a  response  with  a 
large  negative  normal  acceleration  spike  at  t  =  0+. 

It  was  especially  noticeable  in  Flight  Conditions  1,  9, 
and  13.  The  spike  appeared,  with  the  exception  of  Flight 
Condition  9,  in  cases  where  the  g.  was  low,  which  is  rea¬ 
sonable  if  one  considers  the  aircraft  in  such  an  environ¬ 
ment.  In  order  to  achieve  a  quick  rvse  time  at  low 
dynamic  pressure,  the  control  action  of  an  elevator 
deflection  translates  the  aircraft  vertically  downward 
before  rotation  occurs.  With  a  high  g.  the  elevator 
has  greater  pitch  control  effectiveness,  thus  the  aircraft 
does  not  experience  so  great  a  negative  nz.  Flight  Con¬ 
dition  nine  is  affected  more  by  other  factors  and  will  be 
discussed  later. 


Pre-Filter 

To  deal  with  the  negative  normal  acceleration  spike 
problem,  a  pre-filter  was  added  to  shape  the  pilot's 
commanded  input.  The  corresponding  equations  are  presented 
in  Chapter  4.  The  filter  was  not  used  initially  in  the 
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optimization  routine  since  the  fixed-form  SAS  was  chosen 
to  bo  solely  the  feedback  system  in  Figure  7  ,  However, 
the  pre-filter  was  added  after  the  simulation  results 
showed  an  undesirable  spike  in  the  output.  Table  IX 
shows  that  the  change  in  pilot  rating  was  significant 
after  the  addition  of  the  filter.  Since  no  extensive 
SAS  optimization  was  performed  with  the  pre-filter  in 
the  system,  it  is  impossible  to  say  just  what  its  effect 
would  be  on  the  final  closed  loop  response.  Probably  one 
of  two  possibilities  would  result.  Either  the  optimiza¬ 
tion  and  resulting  feedback  gains  would  bring  a  negative 
spike  back  into  the  response,  or  the  filtered  input  would 
be  sufficient  to  yield  a  final  response  ’without  the  high 
frequency  spike.  However,  this  study  was  concerned  only 
with  augmenting'  the  aircraft  with  a  feedback  SAS,  and  the 
filter  was  added  merely  to  check  its  effect  on  the  negative 
spike,  regardless  of  its  effect  on  the  rest  of  the  response. 
It  was  concluded  that  the  original  SAS  design  was  still 
valid  in  its  "optional"  sense. 

The  pre-filter  used  was  a  first  order  lag.  Break 
frequencies  of  both  2  and  5  radians  per  second  were  tried 
for  each  flight  condition. 

The  addition  of  the  filter  brought  the  responses 
closer  to  the  prescribed  envelopes,  but  in  general  the 
negative  spike  was  still  large  enough  to  fall  outside  of 
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the  envelope.  The  effect,  of  the  filter  depended  on  the 
"frequency"  of  the  spike.  A  short  narrow  spike,  such 
as  on  flight  condition  5,  was  effectively  removed  by  the 
filter,  while  a  long,  wide  spike,  such  as  on  flight 
condition  13,  remained  with  the  filter  in  operation.  For 
the  most  part,  the  filter  with  a  break  frequency  of  2 
was  too  restrictive,  lowering  rise  times  below  acceptable 
limits.  Perhaps  a  new  optimization  with  this  filter 
employed  would  yield  better  performance  characteristics. 


High  Free 


Low  Danninc 


The  closed  loop  response  for  flight  condition  9  falls 
outside  the  Design  Challenge  envelopes.  The  reason  for 
the  failure  of  the  Paper  Pilot  in  this  case  seems  to  be 
the  same  problem  as  discussed  in  Chapter  II.  The  Paper 
Pilot  does  not  accurately  reflect  pilot  opinion  in  the 
case  of  an  aircraft  with  a  high  natural  frequency  and  low 
damping.  Table  XII  shows  the  open  loop  short  period 
parameters  for  the  simulated  flight  conditions.  The  final 
column,  "Improvement  ii.  Closed  Loop  Response",  is  a 
judge  of  how  much  the  SAS  improved  any  deficiencies  in  the 
open  loop  response.  Flight  condition  9,  with  a  $  of 
0.109  and  of  5.51,  falls  into  the  region  where  the 

Paper  Pilot  had  previously  failed.  It  must  be  concluded, 
therefore,  that  the  addition  of  remnant,  discussed  in 
Chapter  II,  did  not  completely  solve  the  problem.  Further 


r 


G.GC/ES/73-3 


Table  XII.  Short  Period  Paramo t<- .  s 
and  Closed  Loop  Response  Quality 


support  of  this  conclusion  can  be  found  in  Table  XII. 

Flight  condition  13,  with  an  equally  poor  open  loop  res¬ 
ponse,  yet  a  lower  frequency,  was  improved  significantly 
when  the  SAS  was  added.  The  final  pilot  rating  was  3.97. 
Flight  condition  9,  received  a  better  pilot  rating,  3.27, 
yet  the  closed  loop  response  was  considered  very  poor. 

Flight  condition  5  also  fell  into  this  high  frequency, 
low  damping  region.  While  the  final  closed  loop  response 
was  not  bad,  the  degree  of  improvement  over  the  open  loop 
response  was  slight,  indicating  once  again  the  problem 
with  Paper  Pilot. 

A  possible  cure  for  this  problem  would  be  the  addition 
of  a  normal  acceleration  term  or  pitch  rate  term  (which  .is 
related  to  normal  acceleration)  in  the  cost  function  used 
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tern  be  added  to  the  Paper  Pilot  rating  expression 
to  correct  this  problem. 

Although  a  good  design  procedure,  Paper  Pilot 
proved  to  be  very  costly  in  computer  time.  The  Air  Force 
Flight  Dynamics  Laboratory  has  recently  improved  search 
procedures  to  lever  the  computer  time.  Also,  reducing 
the  state  vector  of  the  entire  system  by  approximations 
significantly  reduces  the  required  time.  In  this  study 
it  was  found  that  ignoring  the  two  second  order  bending 
modes  of  the  Design  Challenge  aircraft  made  little 
change  in  the  SAS  design  yet  greatly  improved  computer 
time . 

Recommendations  for  Future  Study 

In  light  of  various  problems  and  insights  encountered 
during  the  course  of  this  study,  the  following  suggestions 
are  offered  for  future  work  using  Paper  Pilot: 

1.  Incorporate  a  pitch  rate  or  normal  accelera¬ 
tion  term  into  the  pilot  rating  expression 
in  Pitch  Paper  Pilot. 

2.  Optimize  pilot  ratings  with  respect  to  pilot 
lag  as  well  as  lead  and  gain,  to  study  effects 
on  high  frequency,  low  damping  cases. 

3.  Find  optimal  SAS  gains  for  the  JACC  Design 
Challenge  with  a  pre-filter  included  in  the 
system. 

Applying  some  of  these  refinements.  Pitch  Paper 
Pilot  could  prove  to  bo  a  very  useful  tool  in  SAS  design. 


C  r 
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ft/soc,  which  is  a  rather  small  vortical  wind  gust. 

It  is  not  unreasonable  to  find  that  the  aircraft  does 
not. respond  to  this  :-:nall  disturbance.  . 

Comparisons 

The  JACC  Design  Challenge  was  studied  earlier  by 
Sherraru,  Ward  and  Prjol.  Sherrard  used  optimal  control 
techniques,  Ward  used  classical  control  techniques  and 
Paul  used  self-adaptive  control  techniques.  Comparing 
the  simulation  results  of  those  studies  with  the  results 
presented  in  Appendix  C  illustrates  the  success  of  Paper 
Pilot.  Paper  Pilot  optimization  gave  better  compliance 
with  the  Design  Challenge  response  envelopes  than  flight 
conditions  3  and  5  in  Sherrard,  1  and  3  in  Ward,  and  1, 

3  and  13  in  Paul.  In  those  flight  conditions,  Paper  Pilo 
had  better  rise  tirr.es  or  damping  in  the  normal  accelera¬ 
tion  responses. 
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VI .  Conclusions 

Paper  Pilot  offers  a  viable  method  of  SAS  parameter 
optimization  for  SAS  design.  Selecting  SAS  gains  that 
give  the  minimum  pilot  rating  for  a  particular  aircraft 
configuration  results  in  SAS  parameters  that  are  pilot 
rating  optimal.  In  this  study,  this  method  often  resulted 
in  a  closed  loop  response  that  was  as  good  as  or  better 
than  the  response  from  designs  by  classical  or  modern 
control  techniques. 

Initially,  Pitch  Paper  Pilot  was  modified  to  include 
terms  for  pilot  remnant  and  lag  in  the  pilot  model.  This 
was  to  improve  the  accuracy  of  the  predicted  pilot  ratings 
for  flight  conditions  with  a  high  short  period  natural 
frequency  and  low  damping.  This  was  not  the  total  answer, 
as  some  inaccuracy  still  existed  after  the  modifications. 

Paper  Pilot  was  next  used  to  select  SAS  gains  for 
the  JACC  Design  Challenge  aircraft.  An  undesirable  nega¬ 
tive  normal  acceleration  spike  appeared  in  the  analog 
simulation  of  the  closed  loop  aircraft.  Adding  a  pre¬ 
filter  to  shape  the  pilot's  commanded  input  improved  this 
aspect,  but  complete  SAS  parameter  optimization  with  the 
pre-filter  included  was  left  to  future  studies.  The 
problem  with  high  frequency  and  low  damping  conditions 
appeared  again,  and  it  was  suggested  that  a  pitch  rate 
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term  be  added  to  the  Paper  Pilot  rating  expression 
to  correct  this  problem. 

Although  a  good  design  procedure.  Paper  Pilot 

proved  to  be  very  costly  in  computer  time.  The  Air  Force 
Flight  Dynamics  Laboratory  has  recently  improved  search 
procedures  to  lover  the  computer  time.  Also,  reducing 
the  state  vector  of  the  entire  system  by  approximations 
significantly  reduces  the  required  time.  In  this  study 
it  was  found  that  ignoring  the  two  second  order  bending 
inodes  of  the  Design  Challenge  aircraft  made  little 
change  in  the  SAS  design  yet  greatly  improved  computer 
time . 


Recommer.dat i c.is  for  Future  Study 


In  light  of  variou'  problems  and  insights  encountered 
during  the  course  of  this  study,  the  following  suggestions 
are  offered  for  future  work  using  Paper  Pilot: 

1.  Incorporate  a  pitch  rate  or  normal  accelera¬ 
tion  term  into  the  pilot  rating  expression 
in  Pitch  Paper  Pilot. 

2.  Optimise  pilot  ratings  with  respect  to  pilot 
lag  as  well  as  lead  and  gain,  to  study  effects 
on  high  frequency,  low  damping  cases. 

3.  Find  optimal  SAS  gains  for  the  JACC  Design 
Challenge  with  a  pre-filter  included  in  the 
system . 

Applying  some  of  these  refinements.  Pitch  Paper 
Pilot  could  prove  to  be  a  very  useful  tool  in  SAS  design. 
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Appendix  A 

System  Equations  in  State  Variable  Form 

Ai rc ra ft  F ^ a aliens  of  Mot \ on 


The  linearized  longitudinal  short  period  equations 
are  (7:3) 


e 

q 

a 

Since  normal 
in  place  of  anqle 
ation  is  needed. 


=  q 

(22) 

(23) 

-  =."+,j*zi /.  +z°  Sf 

(24) 

acceleration  is  used  as  a  feedback  state 
or  attack,  an  equation  for  normal  acceler- 
This  is 


n  =  ^  (q-i)  (25) 

9 

Substituting  equation  (24)  into  equation  (25)  and  ignoring 
the  wind  gust  terms  yields 


nz  - 


Using  the  fact  that  w  = 

Z."  -  UQZ  )  equations 

««  o. 


q  = 


w  = 
nz  = 


UQ 

—  j.  ) 

b  z«“  - 

Uq a  and  letting  Mw 
(23),  (24)  and  (26) 


My*w  +Mqq+M^e 

z»,w  +u0q+z5  ^ 

°e 


can  be 


(26) 

Zw~  Za  • 

rev.’ritten 


(27) 

(28) 
(29) 


Preceding  page  blank 
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For  this  aircraft,  the  break  frequency  of  the  filter 
was  chcsam  as  3.  This  was  the  value  arrived  at  by 
Sherrard,  end  represented  a  good  choise  of  maximum 
attenuation  cf  body  bonding  effects  with  minimum 
attenuation  of  the  rigid  body  dynamics  (8:28). 

From  the  design  challenge  the  indicated  pitch  rate 
and  indicated  normal  acceleration  measurement  equations 


are 


<3i  =  -A69u2  +qn 


(31) 


n2i  =  nz+g  ^zq+Pn^  +!2n2  V+nzn  (32) 
The  first  order  filter  yields  the  two  expressions 


nZf  -  -w£n2£  +  ^nZ£ 

(33) 

qf  =  -  wf  q  f  +  wfqi 

(34) 

The  first  and  second  body  bending  modes  are 

represented 

by 

2 

Wj  +  (2  5^  )  01  +  (  Uj )  u1  =  f  Se 

(35) 

.2 

^2+  (2J2oj2)  02  +  (  uj^)  u2  =  y  Sc 

(36) 

Letting 

•  d~ 

II 

Vi 

S  =  U2 

(37) 

and  substituting  in  the  values  of  ^  , o  and  y 

yield  the 

body  bending  equations  as  used  here 

r  +  .  06r  +  900  .1  =  -4U0Zj  fe 

(38) 

s  +  s  +  2500  >.2  =  -4UqZ)j  Se 


(39: 
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The  final  SAS  output  is 

js=  K^e  +  Kqqf  +  KnnZf  (40) 

Multiplying  equation  (34)  by  Kq,  equation  (33)  by  Kn, 

and  combining 

Kq4f  +  Kn^*Zf 

=  -tof  (Kqqf  +  Knnz..l  +  uf  ( Kc.qi  +  Knnz.)  (41) 

Now  define  a  new  state. 

/>  4  Kqqf  +  KnnIf  (42) 

Differentiating  this  and  equating  the  result  to  (41) 

ft  =  K„cr-  +  K„n,^ 

r  q  *f  n  H  f 

ft  -  -uf/5  +u.yKqqi  +  wfKnnZi  (43) 

Substituting  equations  (27)  and  (29)  into  (32)  yields 

°Zi  =  'Szww  ~  qzje4 

+  |  (L2Mww+LzMqq+L2M  4  )  (44) 

+  ~  +  +  nZn 

Finally,  after  substituting  equations  (31)  and  (44)  into 
er nation  (43)  ,  the  final  equation  for  ft  results 

ft  “  “  “f  P 

+  Lz-V  q 

+  "fib  <-Z„  *LzM„)w 

+  ['^e  +^"’4  -  4bc'  (IV0"2>]  4 
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Appendix  B 


State  Equations  from  the  Pilot  Model 


As  shown  ir.  Figure  2 ,  the  pilot  modal  responds  to  a 

pitch  error  through  a  lag  and  lead  network  and  pure  time 
delay,  with  remnant  added  at  the  output.  The  output  of 
the  pilot  model,  modified  by  the  system  feedback,  is  the 
command  to  the  actuator. 


Time  Delay.  For  purposes  of  simplicity  of  calcula¬ 
tion  and  modeling,  the  time  delay  is  approximated  by  the 


Pade  expansion,  thu? 


TS  5  !  -s+0.5r) 

e  =  ~y  ~  (s+  0,5/7 


This  can  be  re-written 


Define  a  new  state 


By  substitution 


s'*s  =  7(j-ii 


y,  =  S'*  S 


y,  7  (-  y,  +  2  S') 


Lag-lead  Network.  Now,  from  the  transfer  function 


S_  r,  y p(TrjS>l) 

©f  (TjS  +1) 


letting  0f=©c-0  and  a=9 


5 T  -KpTL9c  =  -$-KDTLq+Kp9c -KD© 


->  I  1NPxL°c  - 
Define  a  new  state 


y2  =Tj  5  '  -KpTLBc 


Preceding  page  blank 
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By  substitution 

y 2  =  -  S  -KpTLq+Kt)9c-Kp9 

Frc-ra  the  equation  to r  y2 

('=  -,34.i’o7’L  o 

j  17 -T7~  c 

Pitch  Command.  From  the  system  diagram,  Figure  7,  the 
commanded  pitch,  ©c  ,  is 

5 


(53) 


(54) 


0,.  = 

c  S  +  1 


(55) 


where  j  is  a  zero-mean  gaussian  white  noise.  Therefore 


6c  =  “^bec+  $ 


(56) 


Actuator  Equation.  Again  from  the  system  diagram, 
the  actuator  equation  is 


f  =  -  1_  (  +  1  c 

VP  -  ZZ — 

rp  0  rp  e 

1e  e 

However,  since,  from  the  system  diagram, 


J-‘JsT'5r  5  ~Ss  +  ST 


157) 


(58) 


and  from  Appendix  A, 


Ss~  Ve  + 


By  substitution 

r  1r  Kr  ,  1  1  f '  1  „  1  /■ 

A  =  -  -zr  j_  -  “y  JG  +  f  y ,  -  —  A  -  —  G  +  -7-  )r 

Remnant.  From  the  system  diagram.  Figure  7 

s  =  — ^ 


(59: 


(oO  J 


'r  s+  w. 


(61) 
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which  can  be  re-written 

4=~"^r+^r  (62) 

System  Equations.  By  proper  substitution  of  the  abov 
results,  the  equations  shown  below  are  formed.  When 
combined  with  the  airframe  and  SAS  equations  from  Appendix 
A,  they  form  a  complete  set  of  system  equations  for  this 
pilot-vehicle  model. 


h 

=  f«-yi^y2+^tec, 

(63) 

Y2  = 

-KPTLq"V"  T?  y2+KP<1*  'T^)S‘ 

(64) 

®c  *  -  -bec+  5 

(56) 

S  = 

-  T«(1+KJ  I'l-  tgTj  i'2 

(65) 

tqtt  KpTLec  +  Te  5r"  Te  P 

ir*  ~Mr  + 


(62) 


c: 


:  :u  e 


X  ry.;  5  o--ii  j:re-::jltv;.=  t  o  t.w'.  ■>r.  the-  t:  lot’s  eoTr.ir-ar.d'c-d 
and  ...  j  r i  •  " t  .  r';  ;  s  '.rc-ii ;  to:  is  a  first  ordc 

It  i  "fr-t  Ui.-i.-d  w_  ‘..h  «  brcaV.  i  resume;/  of  5  rudi- 
second,  t'ton  reduced  to  2  raciar.s  per  second. 


Trie 


filial  plot  is  that  of  the  a  if  era 


u  r 


.cor 


the  specified  vertical  wind  oust  (2.821  ft/soc  rr.'.i 


st  to 

)  .  All 


flight,  conditions  here,  both  open  and  closed  loop,  showed 
zero  response  to  the  gust  on  the  indicated  scale.  There¬ 
fore,  only  one  flight  condition  is  presented. 
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--Filter  Equations 

(66) 
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(/> 

(67) 

•  /  .  ' 

-  +  “p^p 

(68) 

From  Appendix  B, 

1691 

Substituting  equation  (69)  into  equation  (66) ,  and  then 
into  equation  (68) 

“pVi"  ^y2"  9c  +“p4 


_  *^p  +  wp^p 


T^-  Qc  +  ^‘fr  -  UP^P  ^7^ 


The  resulting  new  state  equation  for  <$pis 

^  9‘  +  (71) 

Also  from  Figure  12 

K-\Sr\‘.  <72) 

Substituting  equation  (67)  into  equation  (72)  yields 

(  -  X  f  -  J-S  -  4-/  (73) 

«>e  -  T„  Te  JS  Tp  "e 


«>e  -  Te  "p  Te  S  ie 
Since  from  Appendix  A 


4  =  ^.  +  p 


(59) 

Preceding  page  blank 


we  car.  sutrt:  t  xze  equation  (59)  into  equation  (73)  to 
get  the  state  equation  in  <fe 


X(,  +  KJ>]  t 


which  is 


change  to  the  previous  equation  for  Se 
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